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We discuss signatures of thcrmalisation in heavy ion collisions based on elliptic 
flow. We then propose a new method to analyse elliptic flow, based on multiparticle 
azimuthal correlations. This method allows one to test quantitatively the collective 
behaviour of the interacting system. 



1 Introduction 

A highly debated issue in ultrarelativistic heavy ion coUisions is whether or 
not the colhding system reaches local thermal equilibrium at some stage of 
its evolution. Equilibrium with respect to inelastic collisions strongly con- 
strains the ratios of particle abundances El ( "chemical" equilibrium) , as well 
as phase-space densities, which are obtained by combiniiig informations from 
momentum spectra and two-particle HBT correlations^ on the other hand, 
equilibrium with respect to elastic collisions constrains momentum distribu- 
tions, and implies in particular that they are isotropic in the local rest frame. 
This is the "kinetic" equilibrium, on which we concentrate here. 

Kinetic equilibrium itself has (at least) two facets. One is the equilibration 
between longitudinal and transverse degrees of freedom, i.e., the implication 
that in the local rest frame, longitudinal and transverse momenta are of the 
same order of magnitude. This aspect of thcrmalisation can be discussed from 
first pripjciples at the partonic level,Q and there is now a vast literature on this 
subject.u But experimental signatures deal in fact rather with equilibration 
among the two transverse degrees of freedom. Due to the high Lorentz con- 
traction at ultrarelativistic energies, the typical transverse scale is much larger 
than the longitudinal scale, so that this "transverse equilibrium" is probably 
easier to achieve than "longitudinal-transverse equilibrium" . 

In Sec. 1^, we discuss some experimental signatures of transverse equilib- 
rium, especially elliptic flow, and the relevance of hydrodynamical models in 
this context. In Sec. ^, we present a new method recently developed to obtain 
more reliable measurements of elliptic flow, which are required in order to 
draw definite conclusions on the issue of thcrmalisation. 



2 Elliptic flow: a signature of transverse thermalisation 



2.1 When are hydrodynamical models useful? 

If the colliding system thermalises, its evolution follows the laws of hydrody- 
namics; it is therefore natural to use hydrodynamical models in order to define 
signatures of thermalisation. There are two ingredients in these models: i) an 
equation of state; ii) initial conditions, i.e., energy density, baryon density and 
fluid velocity on a space- like hypersurface (typically, at some initial time). 

Which parameters are under control? The initial time to, at which the 
system thermalises, is hopefully short but to a large extent unknown, so that 
hydrodynamics has little predictive power concerning observables which de- 
pend strongly on to, such as thermal photon productions Similarly, there is 
no serious motivation for studying the very first stages of the collision within 
the framework of hydrodynamics, either within a one-fluid model I3'LI or a mul- 
tifluid model, u and the relevance of signatures based on such parametrisations 
is questionable. Finally, there is no well-defined prescription concerning the 
initial longitudinal fluid velocity and density, foii which a variety of parametri- 
sations exist, either "Landau" □ or "Bjorken" i3 type. 

The situation is much clearer concerning transverse degrees of freedom: 
the transverse collective velocity must be initially zero since each nucleon- 
nucleon collision populates the transverse momentum space randomly; the ini- 
tial density profile in the transverse plane is strongly constrained by observed 
multiplicity distributions JIj Therefore, reliable signatures of thermalisation 
should rather be sought Jpr in observables associated with transverse degrees 
of freedom: px spectra,tj transverse radii t3 and azimuthal anisotropics, in 
particular elliptic flow, on which we concentrate here. 

2.2 Why elliptic flow? 

Elliptic flow is defined as a correlation between the-azimuth (jj of an outgoing 
particle and the azimuth of impact parameter:^ 

^;2 = (e^'(^-*-'), (1) 

where brackets denote a statisticaLaKerage. At ultrarelativistic energies, V2 is 
positive for noncentral collisions .lJIIj In a hydrodynamical picture, it results 
from anisotropic pressure gradients in the transverse plane, due to the almond- 
shaped region of the overlap region between the two nuclei.EJ Microscopically, 
1)2 is created by rescattering among the produced particles, which makes it a 
sensitive probe of final state interactions: if there are none, it vajiishes, while 
other observables such as px spectra may still look "thermal" .113 



Furthermore, predictions of hydrodynamical models for V2 are very stable: 
since V2 is created by transverse pressure gradients, it strongly depends on the 
initial density transverse profile, which is well controlled as discussed above; 
it also depends significantly on the equation of state, on which it may thus 
provide valuable information; on the other hand, V2 depends only weakly on 
arbitrary parameters, such as initial time and longitudinal velocity.liZI Quite 
remarkably, simple hydrodynamical parametrisations are able to reproduce aii 
multaneously the measured spectra, HBT radii and elliptic flow at RHlCJlj 

2.3 Predictions for V2 

The momentum anisotropy V2 calculated in hydro rnodels is roughly equal to 
the anisotropy of the almond-shaped overlap regionO This purely geometri- 
cal effect dominates the centrality dependence of V2 , which decreases linearly 
with the number of particinants. Deviations from this behaviour can ha. used 
to signal a phase transition or a departure from thermal equilibrium Eil The 
latter is expected to occur for the most peripheral collisions, where V2 should 
be smaller than the hydro prediction—This is indeed observed in several spe- 
cific transport models like_IJrQMD c3 (which however predicts a much too 
smaU value of i;2), QGSM,Ej and AMPT.0 With a more systematic studyJ^S 
one could relate the observed centrality dependence of V2 to the degree of 
thermalisation of the system. Experimental results,jiary significantly depend- 
ing on the method used to analyse elliptic flow.E^O We come back to this 
issue in Sec. ^. j— . 

Hydrodynamical calculations E3 were also able to predict the pt depen- 
dence c£^2 for identified hadrons, in remarkable agreement with experimental 
results:Ej V2 is almost linear in pT for pions and significantly smaller for pro- 
tons. IJfmffjycj', these non-trivial features are also reproduced by transport 
models.E3l2j'Eil In addition, the latter predict a saturation at high pT which 
is not seen in hydro, suggesting that many elastic collisions are necessary to 
build the flow at high pT- This_saturation, which is seen in the data, could 
also be related to hard physics.E5l 

3 Analysing elliptic flow with multiparticle correlations 

3. 1 Flow from azimuthal correlations 

Since the reaction plane $ij is unknown, V2 cannot be derived directly from 
(^. It must be inferred from azimuthal correlations. The standard flow 
analysis Ell relies on the key assumption that all azimuthal correlations are due 
to flow, i.e., that angles relative to the reaction plane 0— are statistically 



independent. This allows one to write the two-particle correlation as I 

g2j(0i-02)\ ^ /g2j(0i-0j,.)^2i(0R-02)\ ^ ^,^^,2 



g2«(0R-02)\ („2)2^ (2) 



where brackets denote an average over pairs of particles belonging to the same 
event. One could also use higher order correlations, such as the four-particle 
correlation, which would give 



2i{<j>i+ct>2—ct>3-<t>4) \ _ 



(^^2)^ (3) 



However, such equations are not quite correct, since they neglect various other 
sources of azimuthal correlations ("nmiflow correlations"), which is no longer 
justified at ultrarelativistic energies.cJ 



3.2 Simple illustration of nonflow correlations 

In order to illustrate nonflow correlations, we consider the following example: 
assume that in each event, M/2 pairs of particles are emitted, where both 
particles in a pair have coUinear momenta, but pairs are emitted with random 
orientations. Since pairs are emitted randomly, there is no flow (i;2 = 0), but 
there are azimuthal correlations. In each event, there is a total of M{M — 
l)/2 particle pairs, among which M/2 are correlated, hence the two-particle 
correlation 

1 



2i(0l-02)\ _ 



M- 1 

A similar reasoning yields the four-particle correlation: 

.2.(0,+02-03-04)\ _ 2M{M-2) ^ 



(4) 



(5) 



M{M -1){M -2) (Af - 3) {M - 1) (M - 3) 

Applying Eqs. (^) and (^), or (^) and (^, one would obtain V2 ^ 1/\/M, al- 
though there is no flow: this is the typical order at which nonflow correlations 
spoil the standard flow analysis. 



3.3 Subtraction of nonflow correlations 

The contribution of nonflow correlations can be greatly reduced by combining 
the informations from two- and four-particle correlations. In Eq. (|5|), a 4-uplet 
of particles gives a nonvanishing contribution if it consists of two correlated 
pairs, either (1,3) and (2,4) or (1,4) and (2,3). Subtracting the corresponding 
contributions, one obtains 

/ 2j(0i+02-03-04)\ _ O /„2j(0i-02)\^ _ 1 (a\ 

V / I - (M-l)2(M-3)- ^ ^ 



The l.-h. s. of this equation defines the cumulant of four-particle correlations, 
C2{4}. The r.-h. s. is the contribution of nonflow correlations, of order 
i.e., much smaller than the corresponding contribution to the four-particle 
correlation (|^), of order l/A'P. On the other hand, the contribution of flow 
remains of the same magnitude: from Eqs. (^, ^ and (|^), one obtains 
C2{4} = — (^2)^. The subtraction therefore reduces the relative contribution 
of nonflow effects, and yields a more accurate estimale of the flow. 

This method was recently applied to STAR data.EZi The repilting value of 
V2 are smaller than those obtained with the standard analysis,Ej in particular 
for the most peripheral collisions: this is precisely where nonflow effects are 
expected to give the largest contribution since the multiplicity M is smaller. 
The centrality dependence obtained with this method suggests that departures 
from thermalisation at RHIC may be larger than was previously thought. 

This cumulant expansion can be worked out to arbitrary orders, anc 
allows one to extract the genuine 4-, 6-particle correlations and beyond,^ 
The practical implementation of the method is described in detail elsewhereH^ 
Flow, which is essentially a collective phenomenon, contributes to all orders, 
while the relative contribution of nonflow correlations decreases as the order 
increases. Higher order cumulants therefore provide a unique possibility to 
check quantitatively that azimuthal correlations are indeed of collective origin. 
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